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■ Abstract 

I Recent key observational results on the variation of fine structure constant, the proton to 

. electron mass ratio and the gravitational constant are reviewed. The necessity to substantiate 

the dark sector of cosmology and to test gravity on astrophysical scales is also emphasized. 

!>■ : 

^ ! 1 Introduction 

^ : 

, The question of whether the constants of nature may be dynamical goes back to Dirac [1] who 

I expressed, in his "Large Number hypothesis" , the opinion that very large (or small) dimensionless 

■ universal constants cannot be pure mathematical numbers and must not occur in the basic laws 
. of physics. He stressed that the ratio between the gravitational and electromagnetic forces 
I between a proton and an electron, Gmemp/e^ ~ 10"^" is of the same order as the inverse of the 

age of the universe in atomic units, e^Ho/meC^- He stated that these were not coincidences and 
that these big numbers were not pure constants but reflected the state of our universe. This led 

O I him to postulate that G varies^ as the inverse of the cosmic time. 

' Diracs' hypothesis is indeed not a theory and it was shown later by Jordan [2] (see Ref. [3] for 

■ details) that varying constants can be included in a Lagrangian formulation as new dynamical 
^ ! degree of freedom so that one gets both a dynamical equation of evolution for this degree of 

I freedom and a modification with respect to the equations derived under the hypothesis it is 

j_j ■ constant. Testing for their constancy is a fundamental test of gravitation related to the local 



position invariance. It was also realized [4] that varying constants may be associated to the 
existence of a composition dependent fifth force and thus to a violation of the universality of 
free fall. 

In this review talk, I first emphasize the necessity to test gravity on astrophysical scales. In 
particular, I will discuss what can be learnt from the tests of the constancy of the constants. 
While it is tested with increasing precisions in the laboratory (see Ref. [5] for a recent review), in 
the Solar System and by the study of pulsar timing (see the contribution by G. Esposito-Farese 

^Dirac hypothesis can also be achieved by assuming that e varies as t^^'^. Indeed the choice depends the choice 
of units, either atomic or Planck units. There is however a difference: assuming that only G varies violates the 
strong equivalence principle while assuming a varying e results in a theory violating the Einstein equivalence 
principle. It does not mean we are detecting the variation of a dimensionful constant but simply that either /he 
or Gm^/hc is varying. 
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in this volume), there exist very few tests on astrophysical and cosmological scales [6]. Such 
tests are however required mainly because the cosmic matter budget [7] mostly relies on the 
properties of gravitation, as well as the conclusion that about 96% of the energy density of our 
universe is dark (including dark matter and dark energy). 

I will then review the recent observational and experimental constraints on the variation of 
the fine structure constant, the gravitational constant and the proton to electron mass ration. Let 
us stress once again that we are only able to detect the variation of dimensionless constants (see 
Refs. [8, 9] for discussions and I refer to Refs. [9, 10, 11] for the distinction between fundamental 
parameters and fundamental units, as well as for their role in the formulation of the laws of 
physics) . 

To finish, I discuss some developments concerning the phenomenology of varying constants 
both in the late and early universe. 

2 Gravity on astrophysical scales 
2.1 The dark sector and gravity 

Most cosmological observations now provide compelling evidences that our universe is under- 
going a late time acceleration phase, the interpretation of which is still a matter of debate. In 
particular, the Friedmann equations for a universe filled only with pressureless matter (including 
dark matter) and radiation 

H{z) = HoE{z) with E'^{z) = ^l,{l + zf + ^^,{l + zf + ^%{l + zf (1) 

cannot explain the current data [12]. Various ways to face this fact have been considered but all 
lead to the introduction of new degrees of freedom, referred to as dark energy, in the cosmological 
scenario, either as new components of gravitating matter or as new properties of gravity^. 

In the first approach, it is assumed that there exists new gravitating components, beyond 
the standard model of particle physics, while gravity is supposed to be accurately described 
by general relativity. Many candidates such as a cosmological constant, quintessence [13], K- 
essence [14] have been proposed (see e.g. Ref. [15] for a review). From a cosmological point 
of view, these models are characterized by their equation of state which can be reconstructed 
from the function E{z). Note that most late time observations (such as diameter and luminosity 
distances or the growth factor of cosmic structures) only depend on some combination of this 
function E{z). 

The other route is to allow for a modification of gravity which means that the long range 
force that cannot be screened is assumed not to be described by general relativity. Many such 
models have been considered. For instance, a light scalar field can couple to matter leading to 
scalar-tensor models of quintessence [16, 17]. This scalar field responsible for variation of the 
gravitational constant may also be, depending on its couplings, at the origin of the variation 
of other constants and of a violation of the universality of free fall (see Ref. [3, 18] for details). 
Other possibilities include braneworld models. Higher dimensional models predict that gravity 
should depart from its standard Newton behavior on small scales and up to now this scale is 
constrained to be smaller than 100 — 500 ^m [5]. Among braneworld models, a class has also the 

^Note a third possibility. To infer the existence of this dark energy, we interpret the cosmological observations 
in a given theoretical frame which assumes e.g. symmetries for the spacetime. It may be that the observations 
are not interpreted in the correct frame. 
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feature to allow for deviations from 4-dimensional Einstein gravity on large scales. This is for 
example the case of some multi-brane models [19], multigravity [20], brane induced gravity [21] 
or simulated gravity [22] where gravity is not mediated only by a massless graviton but include 
a tower of massive gravitons. 

The dark sector plays an increasing role in cosmological models. By testing the theory of 
gravity on astrophysical and cosmological scales we will strengthen these conclusions. These 
tests will contribute to substantiate the physics of the dark sector. Dark matter may be more 
complicated than a pure collisionless gas and dark energy may require to go beyond a pure scalar 
field interacting only with gravity. Concerning the dark energy phenomenology, the reconstruc- 
tion of the function E{z) will not be sufficient to distinguish between many models'^. In Ref. [24], 
we proposed a classification of these models and describes some of the specific signatures that 
can discriminate between them. It is recalled in Fig. 1. 

From a theoretical point of view, string theory seems to be the only known promising frame- 
work that can reconcile quantum mechanics and gravity, even though it is not yet fully defined 
beyond the pcrturbative level. One definitive prediction drawn for the low-energy effective 
action is the existence of extra-dimensions and of a scalar field, the dilaton, that couples to 
matter [25] and whose expectation value determines the string coupling constant. It follows 
that the low-energy coupling constants are in fact dynamical quantities. When the dilaton is 
massless (or almost) it leads to 3 effects: (i) a scalar admixture of a scalar component inducing 
deviations from general relativity in gravitational effects, (ii) a variation of the couplings and 
(iii) a violation of the weak equivalence principle. 

From this perspective, testing gravity may also reveal the existence of further gravitational 
fields or of extra-dimensions and it opens an observational window on the low-energy limit of 
string theory and/or on the stabilization of the dilaton and extra-dimensions. 

2.2 What is tested and what should we test? 

Up to now, the observational status concerning the tests of gravity is the following. 

1. On Solar System size, the Newton law as well as the universality of free fall are tested 
with a very good accuracy (see Ref. [5] for a summary of the constraints) . 

2. On galactic scales, there arc a number of astrophysical constraints that a successful mod- 
ification of gravity will have to face (see e.g. Ref. [26]). If the modification of gravity has 
some relevance on galactic scales then it will have to explain the flattening of the rotation 
curves and to account for the dependence of the galaxy rotation curve on the luminosity 
of the galaxy. This dependence is encapsulated in the TuUy-Fischer relation relating the 
luminosity of a spiral galaxy to its asymptotic rotation velocity L oc f^, with a ~ 4. This 
sets severe constraints on theories in which the cross-over scale with standard gravity is 
fixed (see e.g. Ref [27]) and favored theories where this cross-over scale depends on the 
considered galaxy. Roughly, one needs this crossover scale to behave as 

~ VloGM/c (2) 

^Indeed, it is always possible [23] to construct a scalar field potential that will lead to the "observed" E{z) (on 
which most of the observable - diameter and angular distances, growth of cosmic structures,...- depend) so that 
the precise determination of the late evolution history of our universe, even if very constraining, will not allow us 
to determine the true nature of dark energy. 
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Figure 1: Summary of the different classes of models and of the specific tests that can help 
distinguish between them (besides the equation of state and the growth of cosmic structures). 
The classes differ according to the kind of new fields and to the way they couple to the metric 
Qixy and to the standard matter fields. Upper-left class consists of models in which a new kind 
of gravitating matter is introduced, e.g. quintessence, K-essence. In the upper-right class, a 
light field induces a long-range force so that gravity is not described by a spin-2 graviton only. 
This is the case of scalar-tensor theories of gravity. In this class, Einstein equations are modified 
and there may be a variation of the fundamental constants. Note also the possibility for this 
new field to couple differently to baryonic and dark matter. The lower-right class corresponds to 
models in which there may exist massive gravitons, such as in some class of braneworld scenarios. 
These models predict a modification of the Poisson equation on large scales. In the last class 
(lower- left), the distance duality relation may be violated. From Ref. [24]. 

where ^o ~ lO^^m. On the other hand, the compatibility between X-ray and strong lensing 
observations tends to show that the Poisson equation holds up to roughly 2 Mpc [28] . 

3. On cosmological scales, there is at the moment no direct tests of gravity. Indeed the growth 
of cosmological structure can put some constraints but usually the observations entangle 
the properties of the matter and gravity. Both the acceleration of the universe and the 
variation of the constants (if confirmed) may be indications of the break-down of general 
relativity in this regime. 

What should we test then? General relativity is based on Einstein equivalence principle that 
includes three hypothesis: (i) local Lorentz invariance, (ii) local position invariance and (iii) 
universality of free fall. If these hypothesis are valid, it is thought that gravity is a geometric 
property of spacetime (see e.g. Ref. [29]). We may thus aim to test for both the Einstein 
equivalence principle and the field equations that determine the geometric structure created by 
any mass distribution. Some tests in these directions have been proposed in the past years. 
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• Lorentz invariance: it can be tested through the propagation of high energy particles. A 
summary of the constraints can be found in e.g. Ref. [30] 

• Poisson equation: on sub-Hubble scales, the Einstein equations in an expanding spacetime 
reduce to the Poisson equation 

A$ = AirGpa^d (3) 

that relates the gravitational potential to the density contrast. It was recently argued [31] 
that the comparison of galaxy catalogs such as SDSS or 2dF and of weak lensing obser- 
vations give a direct test of this equation and can be extended up to 100 Mpc (see also 
Ref. [32]). 

• Distance duality: as long as photons travel on null geodesies and the geodesic deviation 
equation holds, it can be shown that there is a reciprocity relation, rs = ro{^ + z) between 
the source and observer area distances. Indeed rs cannot be measured so that this relation 
cannot be tested. However, if the number of photons is conserved, it translates to a distance 
duality relation between the luminosity and angular distances, = Da{^ + z)'^ that can 
be tested. This proposition [33] has been tested [24] using X-ray and Sunayev-Zeld'ovich 
observations of galaxy clusters and no departure from the standard expectation has been 
found (see Fig. 2). This can constrain e.g. models involving photon-axion mixing [34]. 

• Growth of cosmic structures: the growth factor of cosmic structures is sensitive to the 
equation of state of the matter that drives the expansion of the universe [35]. It was 

recently proposed to use the skcwness and higher moments of the density field induced 
during the non-linear Newtonian clustering to test for the existence of a long-range Yukawa 
force [36] . Similarly, it was argued that this skewness was sensitive to the coupling of dark 
energy to dark matter [37, 38]. 

• Testing for the constancy of the (non- gravitational) constants is a direct test of the local 
invariance position^ that makes it a test of Einstein equivalence principle and of the strong 
equivalence principle when the gravitational constant is considered. It has been investi- 
gated by many methods and is the subject of the following section. Note that if violated, 
one expects also a violation of the universality of free fall mainly because the self-energy 
is composition dependent and gravitates. 

3 Update on the constraints on the variation of the constants 

Since the claim [41] that the fine structure constant, a, may have been smaller in the past, there 
have been a tremendous increase in the interest of the fundamental constants of nature and to 
whether there are really constant. Various reviews [3, 42] exist and details and more references 
can be found in Ref. [3] 

''The local position invariance implies that the (non-gravitational) laws of physics determined locally take the 
same form at amy spacetime point. This demands that some uniformity is necessary for reasonable predictions to 
be made about distant part of the universe. This can also be called a local predictability assumption [39] . Indeed, 
it (iocs not exclude theories in which the gravitatioual constant varies. Assuming we can determine its law of 
variation, it just implies that local physics is more complex than we might have thought initially. Again, we see 
the interplay between cosmological tests and local physics. 
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Figure 2: rj = Da/Dl{1 + 2;)^ as a function of the redshift for the 18 clusters of the Reese et 
al. [40] catalog. The error bars include the observational error bars as determined by Reese et 
al. [40] and the uncertainties in the cosmological parameters. Outlayers correspond to bimodal 
clusters that cannot be well fitted by a /3-profile. This shows that there is no significant departure 
from the expected distance duality relation. Prom Ref. [24]. 

Wc review the recent developments concerning mainly the constraints on the variation of the 
fine structure constant, as well as the proton to electron mass ration, /x, and the gravitational 
constant, G. 

3.1 Laboratory constraints 

There have been some marked improvements on the constraints on the variation of a in labora- 
tory experiments. These methods are based on the comparison of atomic clocks using different 
types of transitions in different atoms. According to the comparison, one can constrain the 
variation of some combination of fundamental constants. For instance, the hyperfine transition 
frequency of alkali can be approximated by [43] 



where ji is the magnetic moment of the nucleus, jjln the nuclear magneton, the Rydberg 
constant and Fj-d a relativistic function [43] which strongly increases with the atom number Z 
(e.g. dlnFi.ci(Za)/dlna ~ 0.74 for ^""^Cs and 0.30 for ^'^Rb). 

The comparison of hyperfine transitions in ^''Rb and ^^^Cs over a period of about 4 years took 
advantage of this sharp variation to show [44] that dln(i/Rb/i/Cs)/dt = (0.2 it 7.0) x 10~-^^yr~^ 
at lo". Neglecting possible changes in the amplitude of the weak and strong interactions and 
thus in the nuclear magnetic moments, it translates to 



Another experiment [45] comparing an electric quadrupole transition in ^^^Hg+ to the 
ground-state hyperfine splitting of ^'"''^Cs over a 3 years period showed that |dln(i/Hg/^'Cs)/dt| < 
7.0 X 10~^^ yr~^. This constrains the time variation of 5Cs("^e/"ip)ct®" so that 



UN rrip 



(4) 



a/a = (-0.4 ± 16) x 10"^^ yr 



-1 



(5) 



a/a\ < 1.2 X lO-^Vr 
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Figure 3: In a theory in which gravity switches from a standard four dimensional gravity to a 
five dimensional gravity above a crossover scale of Vg = ^Qh~^ Mpc, there axe different cosmo- 
logical implications concerning the growth of cosmological perturbations. Since gravity becomes 
weaker on large scales, fluctuations stop growing [left panel]. It implies [right panel] that the 
density contrast power spectrum (thick line) differs from the standard one (thin line) but, more 
important, from the gravitational potential power spectrum (dash line). Prom Ref. [31]. 

if both the gyromagnetic factor qq^ and me/nip are assumed constant. 

The comparison [46] of the absolute IS" — 25" transition in atomic hydrogen to the ground 
state of cesium combined with the results of Refs. [44, 45] yields the two independent constraints 

a/a = (-0.9 ±2.9) x 10-^Vr"\ dln(fRb/i^Cs)/dt = (-0.5 ±1.7) x IQ-^^yr-^ (7) 

The comparison [47] of optical transitions in ^^^Yb"^ to a cesium atomic clock at two times 
separated by 2.8 years has shown that dln(z/Yb/^^Cs)/dt = (—1-2 ± 4.4) x 10~^^yr~-^ which 
translates to 

a/a = (-0.3 ± 2.0) x 10"^^ yr'^ (8) 

These methods allow to set very sharp local constraints and, as illustrated by the results of 
Ref. [46], they can be combined to set independent constraints on various constants. 

3.2 Geochemical constraints 

Sharp constraints on the time variation of the fine structure constant can be also obtained from 
the Oklo phenomenon and from the study of the lifetimes of long-lived nuclei. 

The Oklo phenomenon is a natural nuclear reactor that operated during 200,000 years ap- 
proximativcly two billion years ago, that is at a redshift z ~ 0.14. The isotopic abundances 
of the yields give access to informations about the nuclear rates at that time. One of the key 
quantity measured is the ratio gl^Sm/gl^Sm of two light isotopes of samarium which are not 
fission products. This ratio of order of 0.9 in normal samarium, is about 0.02 in Oklo ores. This 
low value is interpreted by the depletion of g^^Sm by thermal neutrons to which it was exposed 
while the reactor was active. The capture cross section of thermal neutron by ^^^Sm 

^^'^Sm + n^^™ Sm + 7 (9) 
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has a resonant energy Er — 0.0973 cV, which is a consequence of a near cancellation between 
electromagnetic and strong interactions [48]. A detailed analysis and model of the samarium 
nuclei lead [49], assuming that the variation of Ej. is due only to the a dependence of the 
electromagnetic energy, to the constraint 

Aa/a = (0.15 ± 1.05) x 10"'^ (10) 

at 2a. In particular, the accuracy of the method can be understood by comparing the resonant 
energy Er ~ 0.1 eV to its sensitivity to a variation of a, dEr/dlna ~ —1 Mev so that variation 
smaller than Aa/a ~ O.leV/lMeV ~ 10~^ are expected. 

It was later pointed out [50] that there may be two ranges of solutions compatible with the 
Oklo data 

Aa/a = (-0.8 ± 1.0) X 10"^ Aa/a = (8.8 ± 0.7) x 10"^ (11) 

the second branches being disfavored by the analysis of the isotopic ratio of gadolinium. 

Recently, the assumption that the low energy neutron spectrum is well described by a 
Maxwell-Boltzmann distribution was investigated [51]. The effect of a variation of the strange 
quark mass was studied [52] to show that ttt-s/Aqcd has varied by less than 1.2 x 10"^'^ while 
assuming all fundamental couplings to vary independently led [53] to the more stringent limit 
Aa/a < (1 - 5) X 10"^°. 

Radioactive decay lifetimes can also be used once the a-dependence of the decay rate is 
known. For instance, the lifetime of a /3- decay nuclei scales as 

A ~ A{AE)PGla' (12) 

where AE is the decay energy and s the sensitivity. Many nuclei were used but the sharpest 
constraint was obtained from the /3-decay of rhenium to osmium by electron emission 

^fRe-^^fOs + Pe + e-, (13) 

first considered by Peebles and Dickc [54]. Interestingly, due its low decay energy -about 2.5 
keV- the sensitivity of rhenium is s ~ —18, 000 so that a variation of a of order 10~^% induces 
a variation decay energy of order of the keV. 

The analysis of new meteorite and laboratory led to [55] 

Aa/a = (8 ± 16) x 10"'^ (14) 

over the last 4.5 Gyr, which corresponds to z ~ 0.45. 

There is a caveat to this method that is not so direct: the ratio Re/Os is measured in iron 
meteorite the age of which is not determined directly. Models of formation of the solar system 
tend to show that iron meteorites and angrite meteorites form within the same 5 million years. 
The age of the latter can be estimated from the ^'^^Pb-^'^^Pb method which gives 4.558 billion 
years. Besides, this constraint holds on the averaged value of a over the 4.5 past billion years. 

3.3 Cosmological constraints 

On cosmological scales, the observation of the cosmic microwave background (CMB) anisotropies 
and of the abundances of the light elements produced during the big-bang nucleosynthesis (BBN) 
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allow to set constraints of order 10 ^ on the variation of a. 

Changing the fine structure constant modifies the strength of the electromagnetic interaction 
and thus its only effect on CMB anisotropics arises from the change in the differential optical 
depth of photons due to the Thomson scattering, r = x^n^ca^^ which enters in the collision 
term of the Boltzmann equation describing the evolution of the photon distribution function 
and where is the ionization fraction (i.e. the number density of free electrons with respect to 
their total number density n^. The first dependence of the optical depth on the fine structure 
constant arises from the Thomson scattering cross-section given by o"t = (87r/3)(?i^/m^c^)Q!^. 
The second, and more subtle dependence, comes from the ionization fraction. A variation of the 
fine structure constant can thus be thought of as considering a delayed recombination model. 

Early works [56, 58, 57] based on BOOMERanG and MAXIMA data tend to show that 
the fit to CMB data are improved by allowing Aa ^ while Landau et al. [59] concluded 
from these data imply, assuming spatially fiat models with adiabatic primordial fiuctuations, 
that -0.14 < Aa/a < 0.03 at 2a level. The recent analysis [60] of the WMAP data (see the 
contribution by G. Rocha for details) gave the 

Aa/a = (-1.5 ± 3.5) x 10"^ (15) 

at z ~ 10^. Note that he variation of the gravitational constant can also have similar effects on 
the CMB [61]. In conclusion, constraints of order 1% on the variation of a can be obtained from 
the CMB only if the cosmological parameters are independently known. 

BBN theory predicts the production of the light elements in the early universe the abundances 
of which rely on a fine balance between the expansion of the universe and the weak interaction 
rates which controls the neutron to proton ratio at the onset of BBN. Basically, the abundance 
of helium-4 is given by 

y = 2 (n/p)fCxp(-tN/r) 
^ l + (n/p)fexp(-tN/r) ^ ' 

where {n/p)f = exp {—Q/kTf) is the neutron to proton ratio at the freeze-out time determined by 
G'p{kTi)^ = s/GN (kT^)^, N being the number of relativistic degrees of freedom; Q = — rup, 
T is the neutron lifetime, Gp the Fermi constant and the time after which the photon density 
becomes low enough for the photo-dissociation to be negligible. As a conclusion, the predictions 
of BBN involve a large number of fundamental constants. 

A change in a affects directly Q and was modelled [62, 63] as Q/Aqcd = aa+bv / Aqcb where 
V determines the weak scale and a and b are two numbers. Roughly, this implies that AY/Y ~ 
—AQ/Q ~ 0.6Ao;/a. On this basis, one can set the constraint lAa/aj < 5 x 10~^ [62, 63], 
confirmed by a recent analysis [64] which gives |Aa/a| < 6 x 10~^. This does not take the 
effect of the fine structure constant in the Coulomb barriers and in the cross-sections, which 
was investigated in Ref. [65]. Recently, the effect of seven parameters {G,a,v,me,T,Q, B^), 
later related to the six constants {G,a,v,me,mu,m(i) was taken into account and led to the 
constraint [66] 

Aa/a = (6±4) X lO-''. (17) 

The effect of the strange quark mass was also investigated [52] and it was claimed that rris/ Aqqd 
has varied by less than 6 x 10~^ since BBN. 
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3.4 Astrophysical constraints 

Most of the excitement over the possibihty of the fine structure constant variation arises from 
the observation of distant quasar absorption systems. The method is based on the comparison 
of absorption spectra to laboratory spectra. 

Initially, the method was based on alkali doublets, the splitting of which gives access to the 
fine structure constant, Aiz/P oc o? . The analysis of Si IV on 21 systems gave [69] Aa/a = 
(—0.5 lb 1.3) X 10~^ for 2 < z < 3. The most recent constraint [70] has been obtained from the 
analysis of SilV in 15 systems of the VLT/UVES sample, improving the former constraint by a 
factor 3, 

Aa/a = (0.15 ± 0.43) x 10"^ 1.59 <z< 2.92. (18) 

It is to be noted that none of the analysis based on the alkali doublet method exhibit a hint of 

variation of a. 

The many multiplet (MM) method proposed by Webb et al. [41, 68] was aimed at increasing 
the precision of the AD method by correlating several transition lines from various species in 
order to reach a sensitivity of order 10~^. In particular, one can compare line shifts of element 
which are sensitive to variation in a with those that are not. At low rcdshift (z < 1.8), the results 
lie mainly on the comparison of Fe to Mg while at higher redshift they lie on the comparison of 
Fe and Si. The latest analysis [67] of the Keck/Hires data (see Fig. 4) based on 128 systems in 
the range 0.5 < 2; < 3 points toward a lower value of a in the past 

Aa/a = (-0.54 ±0.12) X 10-^ 0.5 < z < 3 (19) 

as previous analysis did [41, 68]. A detailed budget of the errors was done and a possible 
systematic effects was looked at but none was exhibited. Note also that new synthetic atomic 
spectra were produced [71]. 

Recent observations from the VLT/UVES Tising the same MM method have not been able 
to duplicate this result [72, 73] (see Fig. 4). The analysis by Chand et al. [72] is mainly based on 
the analysis of Fe and Mg in 23 systems toward 18 QSOs in the range 0.4 < z < 2.3 because they 
apply some selection criteria to the lines. In particular (1) they kept only species with similar 
ionization potentials (Mgll, Fell, Sill and AlII) so that they are most likely to originate from 
similar regions in the cloud, (2) absorption lines that are contaminated by atmospheric lines were 
rejected, (3) they put a threshold on the column density so that all Fell multiplets are detected 
at 5a, (4) they checked that the anchor axe not saturated (Mg I and II are fairly insensitive 
to variation of a) and (5) they excluded strongly saturated systems with large velocity spread. 
They concluded [72, 73] that 

Aa/a = (-0.6 ± 0.6) x 10"^ 0.4 < z < 2.3 (20) 

at la. The analysis of a single quasar [74] also gave 

Aa/a = (-0.1 ± 1.7) X 10-^ z = 1.15 (21) 

mainly from the analysis of Fe lines. These results were further confirmed by the analysis of 
Fe II lines in an absorption system at z = 1.839 toward the quasar QllOl-264 [75] 

Aa/a = (2.4 ±3.8) X 10~^ 2 = 1.839. (22) 
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The combined Fe II sample also gives ls.a/a = (0.4 it 1.5) x 10 ^ for two systems at z = 1.15 
and 1.839. 

Both results [67, 72] are sensitive to the isotopic abundances of magnesium and assume a 
solar ratio ^^Mg:^^Mg:^^Mg= 79:10:11. It is however commonly assumed that heavy magnesium 
isotopes are absent in low mctallicity environments such as absorption clouds. Assuming a ratio 
1:0:0 would have led to a more sig nificant detection of Act/ a = (-0.98 ± 0.13) x 10"^ and 
Aa/a = (-0.36 ± 0.06) x 10"^ respectively for the Keck/Hires [67] and VLT/UVES [72] data. 
This has led to a new interpretation [76] of the MM results in which the variation of a is explained 
by an early nucleosynthesis of ^^Mg and ^^Mg. Interestingly a ratio [72] (^^Mg+^^Mg)/^^Mg = 
0.62 ± 0.05 and 0.30 ± 0.01 can explain respectively the Keck/Hires [67] and VLT/UVES [72] 
data. A model of nucleosynthesis in which ^^Mg and ^^Mg are produced by intermediate mass 
(4 — 6M0) in their asymptotic giant branch was proposed. They can reach a temperature larger 
than 7 x lO^K so that proton capture processes in the Mg-Al cycle are effective enough. This 
hypothesis can be tested by looking at other heavy elements produced by these intermediate 
mass stars. 

To finish, let us mention the use of O III emission lines of quasars and galaxies. The 
analysis [77] of 42 quasars from SDSS early data release and of 165 quasars of SDSS data release 
1 led respectively to the two constraints 

Aa/a = (0.7 ± 1.4) x 10"^, Aa/a = (1.2 ± 0.7) x lO^'' (23) 

in the range 0.16 < z < 0.8. The use of the OH microwave transition was also proposed [78] 
and the preliminary analysis of the quasar PKS1412+135 led to [79] 

Aa/a = (0.51 ± 1.26) x 10"^, z = 0.2467. (24) 

These methods can be extended to higher redshifts {z ^ 3 — 5) and have different systematics 
compared to the MM and AD methods which makes them complementary. Also, the analysis of 
OH, combined with HCO"*" lines gave the simultaneous bounds [80] Aa/a = (— 0.38±2.2) x 10~^, 
An/n = (-0.27 ± 1.6) X 10-3 and Agp/gp = (-0.77 ± 4.2) x lO'^ at z=0.68. 

3.5 Other constants 

3.5.1 Proton to electron mass ratio 

The observation of vibro-rotational transitions of H2 in damped Lyman-a systems allows to 
constrain the variation of the proton to electron mass ratio, ji = m-p/rrie- The spectral lines of 
the quasar are usually translated into a reduced redshift, defined as C,i = {zi — z) / {\ + z) where 
z is the averaged redshift of the absorption system. It can be shown that C,i = {Ajj,/ fx)Ki where 
the Ki are sensitivity coefficients. 

Various constraints have been obtained since 1975 all showing no hint of variation. In 
particular, the analysis of 83 absorption lines [81] gave the limit 

Aji/ij, = (-7.5 ± 9.5) X 10"^ (25) 

at a 2a level. More recently, the analysis of vibro-rotational lines of molecular hydrogen for the 
two quasars Q1232+082 {z = 2.3377) and Q0347-382 {z = 3.0249) tend to show [82] that 

Afi/n = (-5.7 ± 3.8) X 10-^ (26) 
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Figure 4: The analysis of quasar spectra on the variation of the fine structure constant from 
(left) the Keck/Hires data [67] and (right) the VLT/UVES data [72]. The left plot depicts 
their previous low redshift (open circles), previous high redshift (open triangles) and new (open 
squares) samples. The raw results with la error bar are shown on the top panel while the 
middle panel shows the results with an arbitrary bining and the bottom panel combines the 
three samples. The right plot present the results from the VLT/UVES data [72]. The dash lines 
represent the weighted mean and la range from the analysis of Ref. [67] while the dash region 
marks the weighted mean and its 3a error bars. 

at 1.5a. It has recently been revised [84] to 

A/i//i = (2.97 ± 0.74) X 10~^ (27) 

It was shown that the limiting factor of the analysis was the precision of the determination of 
the spectra in the laboratory. New data for the transition wavelengths of H2 Lyman and Werner 
bands have been obtained with an accuracy of 5 x 10"*^ [83]. The reanalysis of the published 
spectra then led to A^//i = (—0.5 it 3.0) x 10^^ at 2a level, confirming, once again, that the 
determination of the laboratory spectra is the key of the debate on a possible variation of /x. 



3.5.2 Gravitational constant 

Few new works concern the gravitational constant. Let us note a new analysis of the BBN [85] 
that tends to show that 

AG/G = 0.01^[!:?^ (28) 

at 68% C.L. It was also shown [64] that the variation of G is correlated to the extra- number of 
relativistic degrees of freedom through AG/G = 78N/A2>. It follows that the combined analysis 
of new '^He and WMAP data implies 

-0.10 < AG/G < 0.13 (29) 
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A recent analysis of the secular variation of the period of nonradial pulsations of the white dwarf 
G117-B15A shows [86] that < GjG < 4.0 x lO^^^ yr^^ at 2a, which is of the same order of 
magnitude of previous independent bounds (see also Ref. [87]). 

More important, in the Solar System, the measurement of the frequency shift of radio photon 
to and from the Cassini spacecraft improved the constraint on the post-Newtonian parameter 
7 to [88] 7 — 1 = (2.1 zb 2.3) x 10~^. This can be translated, in a model dependent way, into 
a constraint on the time variation of G. For instance in a Brans-Dicke theory in a matter 
dominated universe, it implies \G/G\ < lO^^^yr"-*^. 

4 Theoretical motivations and modelling 
4.1 From string theory to phenomenology 

Most higher dimensional theories, such as Kaluza-Klein and string theories, imply that dimen- 
sionless constants are dynamical [25, 89]. For example, in type I superstring, the 10-dimensional 
dilaton couples differently to the gravitational and Yang-Mills sectors because the graviton is 
an excitation of closed strings while the Yang-Mills fields are excitations of open strings. For 
small value of the volume of the extra-dimensions, a T-duality makes the theory equivalent 
to a 10-dimensional theory with Yang-Mills fields localized on a D3 brane. When compact- 
ified on an orbifold, the gauge fields couple to fields Mj living only at these orbifold points 
with coupling q which are not universal. Typically, one gets that = e^^^YgMf while 
S'yM ~ e~'^^VeMj + CiMi. Loop corrections have also been studied in heterotic theory by in- 
cluding Kaluza-Klein excitations [90]. In the limit where the volume is large compared to the 
mass scale, gyM — ©"^^Vg-^fl- ~ baiRMn)'^ /2 + . . .. Again, they arc not universal. It follows 
that the 4-dimensional effective couplings depend on the version of the string theory, on the 
compactification scheme and on the dilaton. 

Interestingly, while the constraints presented above assume that only a was varying, it is 
to be expected that if it varies then all other constants also do. In the context of unified 
theories, it is possible to derive relations between the variations of various constants, which can 
be used to derive sharper observational constraints. The dominant effects [97] arise from the 
variation of the QCD scale, Aqcd and the weak scales v. In particular it was argued [91, 92] 
that A Aqcd /Aqcd ^ 30Aa/a and Av/v ~ 80Aa/a. 

Many phenomenological models starting from the investigation by Bekenstein [93] have been 
developed [94, 95] as well as some braneworld models with varying constant were constructed [96]. 
Damour and Polyakov [97] proposed to capture the features of loop corrections by modelling 
them as a genus expansion, that is as a scries in the string coupling constant. The low energy 
action involves various couplings of the effective four-dimensional dilaton to the different matter 
fields. It follows that generically there appear couplings to matter fields via their Yang-Mills 
couplings and masses as as a potential for the scalar field. This construction includes others 
such as the Bekenstein model [93]. Note that these models are required in order to compare 
constraints at various time. 

This phenomenology is related to the one of quintessence and the light field responsible for 
the time variation of the constants may also be the cosmon [18]. In that sense, the variation of 
the constant may shed some light on the physics on dark energy. Initially models of scalar-tensor 
quintessence [16, 17] were considered. Some were using the Damour-Nordtvedt attraction mech- 
anism toward general relativity to pass the Solar System constraints. The Damour-Polyakov 
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model [97] was generalized [98] to a runaway potential so that the light dilaton accounts for the 
variation of the constants, the acceleration of the universe and is at the origin of a violation of 
the universality of free fall. 

4.2 Two Dangers 

The construction of phenomcnological models that tend to explain the small drift of the constant 
in the late universe by introducing a slow-rolling scalar field have two dangers to avoid. 

First, such a light field will obey a Klein-Gordon like equation, ^ -|- 3H<p = —m?(f) -|- . . .. In 
order for this field not to oscillate but still be evolving, its mass needs to be very small, typically 
Ho ~ 10-^^ eV. The question arises of the mechanism that protects it from radiative 
corrections, a problem common with most quintessence models. Various solutions have been 
proposed among which the possibility for this field to be a pseudo-Goldstone boson [99] or to 
identify this light field with shape modulus [100]. 

A second danger lies in the violation of the universality of free fall due to the composition 
dependence of the self energy and of the masses. This was illustrated in Bekenstein original 
construction [93] and was further studied with linear [102] and quadratic [101] couplings. In 
the case of a light dilaton, it was shown that if it were to remain massless then it would in- 
duce a violation of the universality of free fall seven order of magnitudes larger than the actual 
bounds [97]. To avoid such a catastrophe, it has either to suddenly take a mass larger than a few 
meV (so that gravity will be compatible with Einstein gravity above a millimeter) or decouple 
from matter [97]. This latter mechanism is analogous to the original Damour-Nordtvedt attrac- 
tion mechanism [103]. Both mechanisms have different implications concerning the variation of 
the coupling constants. A fixed point [104] or the recently proposed chameleon mechanism [105] 
also claimed to bypass this problem. A consequence is that the improvement of the tests of the 
universality of free fall by 2 or 3 orders of magnitude may give some surprises. 

4.3 Varying constants in the early universe 

Inflation universally produces classical almost scale free Gaussian inhomogeneities of any light 
scalars. Assuming the coupling constants at the time of inflation depend on some light mod- 
uli fields, it was shown [106] that modulated cosmological fluctuations are produced during 
(p)reheating. This idea was extended to hybrid inflation [107] where the bifurcation value of 
the inflaton is modulated by the spatial inhomogeneities of the couplings. As a result, the 
symmetry breaking after inflation occurs not simultaneously in space but with the time laps 
in different Hubble patches inherited from the long-wavelength moduli inhomogeneities. In this 
model, the consistency relation of inflation is modified. These light field can also be at the origin 
of non-Gaussianity [108]. 

If couplings depend on the value of some light fields, then they have most probably developed 
super-Hubble correlations of typical amplitude 10~^, simply because the quantum fluctuations 
of any light field are amplified during inflation. It follows that one expects spatial variation 
with these correlation (see e.g. Ref. [98]). They may have some observational effects [109], in 
particular on the CMB [110] polarization and Gaussianity. 
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5 Conclusions 



In conclusion, testing for the variation of fundamental constants is a test of fundamental physics, 
and in particular of general relativity. It completes other tests that can be applied on cosmo- 
logical scales. Such tests are needed to substantiate the physics of the dark sector that plays an 
increasing role in cosmology. 

While observational constraints become sharper, the debate on a is not over yet but recent 
observations have not been able to reproduce the detection of a lower o. in the past [41, 68, 67]. 
More important, there now exists a physical model, requiring no new physics, to interpret the 
supposed variation of a as an enhancement of heavy isotopes of magnesium. We can now hope 
that the debate will be settled in the coming years. 

The future will offer new constraints, particularly with atomic clocks in space (ACES), new 
tests of the universality of free fall (recent launch of ^SCOPE) as well as new methods, as 
illustrated by the recent activity. These developments will probably shed some light on a pos- 
sible scalar field acting in the late universe or on some new structures such as higher dimensions. 
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